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Abstract. - A new vibrational level of the Hj molecular ion with binding energy of 1.09 x 
10 _9 a.u.~ 30 neV below the first dissociation limit is predicted, using highly accurate numerical 
nonrelativistic quantum calculations, which go beyond the Born-Oppenheimer approximation. 
It is the first excited vibrational level v = 1 of the 2pa u electronic state, antisymmetric with 
respect to the exchange of the two protons, with orbital angular momentum L — 0. It manifests 
itself as a huge p-H scattering length of a = 750 ± 5 Bohr radii. 



The HJ" molecular ion was one of the first non-trivial quantum mechanical systems treated 
in the early days of Quantum Mechanics 1 . Since, its study has been constantly pursued and 
the bound state calculations reach nowadays a very high degree of accuracy [2I21SE1E1IZI • 

The Born-Oppcnhcimcr approximation provides us with a simple, although approximate, 
description of the system. Two electronic energy curves are correlated with the first dissocia- 
tion limit: the ls<7 g curve with a rather deep energy well supporting twenty L = vibrational 
states (L is the orbital angluar momentum) and the 2po~ u which is mainly repulsive. The later 
one presents however at large distances a weakly attractive potential whose depth is about 
6 x 10~ 5 a.u. A single L = level (v = 0) with binding energy 1.56625 x 10 -5 a.u. has been 
found in this well up to now. 

We will show in this letter that an excited L = level (v = 1) of the 2per„ state of H^ 
actually exists. Its binding energy turns out to be extremely small - 1.09 x 10~ 9 a.u., i.e. 
30 neV or 0.00024 cm -1 - and the corresponding wave function has a spatial extension of 
several hundred Bohr radii, a noticeable fraction of a micrometer. This state manifests itself 
in a huge p-H scattering length a = 750 a.u. which dominates the low energy scattering 
cross section of proton by atomic hydrogen. The Hj formation rate, as well as the subsequent 
abundance of H2 molecules, can be substantially influenced by this resonant p-H cross section. 
The possible existence of an excited vibrational L = level was questioned in [3] 

Calculations have been performed in the framework of non relativistic dynamics and in two 
different ways. On the one hand, by solving the Faddeev equations 8 , what provides the low 
energy p-H scattering parameters and - by mean of the extended effective range expansion 9_ 
- allows to determine the binding energies of near-threshold H^~ states. On the other hand, 
by using adhoc bound state methods [7j which provide the best existing binding energies of 
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the ITj~ system. Despite their radical differences, in their formal as well as in their practical 
implementations, these two approaches well agree in their bound state predictions. 

We summarize in what follows the scattering and bound state results obtained by these two 
methods. As our calculations do not involve Born-Oppenheimer approximation, we use only 
exact quantum numbers in non-relativistic dynamics, where spin-orbit and spin-spin couplings 
are neglected (their effects are discussed below) i.e.: the total orbital angular momentum L, 
the parity II and the symmetry property of the spatial wave function with respect to the 
two-protons exchange; for the latter symmetry, the symmetric states have a total proton spin 
S pp = (singlet) and the antisymmetric ones have S pp = 1 (triplet). The connection between 
exact and approximate quantum numbers is discussed in for L = 0, the lscr s levels are 
singlet, even parity states, while the 2pa u levels are triplet, even parity states. 

The 3-body (p,p,e~) Faddeev calculations are performed in configuration space |SJ. Three 
sets of Jacobi coordinates, corresponding to the different asymptotic states, are involved and 
defined by 



2mpm 1 2m a (mp + m 7 ) rnpTp + m 7 r 7 
-(!•£- r 7 ) y<x = \ — — — — (r tt — ; ) 



y mp + m 7 y m a + mp + m 7 mp 

where (af3j) denote cyclic permutations of (123), m a the particle masses and we identify l=p, 
2=p, 3=e~. The standard Faddeev equations read 

(E-H - V a )^ a = V a J2 */3, (1) 

where Hq is the 3-particle free hamiltonian and V a the 2-body Coulomb potential for the 
interacting ((3"f) pair. 

Equations Q provide satisfactory solutions for bound states but are not suitable for scat- 
tering Coulomb problems. The reason is that their right hand side does not decrease fast 
enough to ensure the decoupling of the Faddeev amplitudes in the asymptotic region and 
to allow unambiguous implementation of boundary conditions. In order to circumvent this 
problem, Merkuriev ^U] proposed to split the Coulomb potential V into two parts by means 
of some arbitrary cut-off function \ 

7(x) - ^(x,y)+^(x,y) 
V s (x,y) = V(x) X (x,y) 
^(x,y) - y(x)[l- X (x, y )] 

and to keep in the right hand side of equation JJJ only the short range V s contribution. One 
is then left with a system of equivalent equations 

(E-H -W- V^ a = V* ]T Vp, (2) 
in which W is a 3-body potential containing the long range parts: 

w = v*+v£ + v* 

The systems of equations Q and J5J are strictly equivalent to the Schrodinger equation and 
realize two different partitions of the total three-body wave function 



* = + * 2 + ^3 = *1 + *2 + *3 
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This approach was found to be very efficient and accurate in calculating the low energy e + Ps 
and e + H cross sections 

Equations J5J) were solved by expanding amplitudes in the bipolar harmonics basis 

* Q (x Q ,y Q )= ^2 y a ,i a {x a ,y<x) B^ M (x a ,y a ) (3) 

and their reduced components (f a ,i a in the basis of two-dimensional splines. Some care has to 
be taken in extracting the scattering observables, specially at zero energy, from the asymptotic 
solution at finite distance. The long range polarization force - generated by the equations 
themselves - makes the convergence of the observables as a function of the p-H distance 
very slow and requires an appropriate extrapolation procedure. The first results of these 
calculations can be found in [8] and a detailed explanation of the method will be given in a 
forthcoming publication [T2] . 

For the lser s , L — state, we obtain the scattering length a = —29.3 a.u. The zero energy 
Faddeev amplitude has 20 nodes in yi-coordinate, indicating the existence of 20 vibrational 
energy levels for H^. 

In the 2per u , L — state, our calculations gives a = 750 ± 5 a.u. The sign and the nodal 
structure of the Faddeev amplitudes indicates that such a big value is due to the existence of 
a first excited state with extremely small binding energy. By calculating the p-H phase shifts 
at small energies and using the effective range theory we are able to determine its binding 
energy. In presence of long range p-H polarization potential, the effective range expansion has 
the form |5] 

k cot S = h aik + a 2 k 2 log k + a 3 k 2 + o(k 2 ) (4) 

where k is the wave vector and ai arc coefficients depending on the interaction. It turns out 
that, in the presence of a weakly bound state with imaginary momentum k — iko, the <xi and 
a 2 terms in (0} become negligible and the expansion recovers the standard form |9"lll3|: 

k cot 6 = a + (3k 2 (5) 

with a — —ko + fik 2 ,. Coefficients a and (3 are obtained by fitting the phase-shift and they 
determine the ko value. Using B = we found |S] by this procedure a bound state at 
B = 1. 13 ± 0.04 x 10~ 9 a.u. below the first p+H dissociation threshold. To our knowledge, 
this is the weakest bond ever predicted, three times smaller than the 4 He atomic dimer |14j . 

The S-wave p-H cross section for the triplet (2pa u ) state is displayed in fig. ^ (dotted 
values). The singlet (lsa g ) contribution is negligible in the zero energy region. It is interesting 
to compare the 3-body calculations with those (solid line) provided by the simple Landau's 
two-body potential V(r) — 2r exp (— r — 1) — 9 /4r 4 ^H] . This model - based on Pauli repulsion 
between protons overbalanced at r ~10 a.u. by attractive polarization forces - gives quite 
a good result for the ground state (B = 1 .4531 x I0~ 5 a.u. instead of the exact 1.56625 x 
10~ 5 value [Jj) but differs strongly in the predictions of the first excited one (B — 4.0566 x 
10~ n a.u.). In the zero energy scattering, both calculations differ also by more than one order 
of magnitude while at energies E ^10~ 8 they are already in quite a good agreement. 

The positive sign of the p-H scattering length and the existence of one node in the Faddeev 
amplitude (see fig.0) unambiguously shows the existence of an excited bound state. After this 
calculation was done, we decided to use a direct method in order to obtain a more accurate 
value of the binding energy as well as a direct computation of its wavefunction. Because 
the three particles are bound, the wavefunction must decrease exponentially if any of the 
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Fig. 1 - Low energy cross section for the scattering of a proton by an hydrogen atom in its ground 
state, for the pp triplet state, compared to the results using the Landau potential for the interaction. 



three inter-particles distances goes to infinity. The idea is thus to expand the full 3-body 
wavefunction on a convenient discrete basis set and to diagonalize the 3-body Hamiltonian 
in this basis set. For highly accurate calculations of very weakly bound states, the basis set 
must be chosen carefully. The first step is to isolate the angular dependance of the 3body 
wavefunction, which is straightforward for L = states. One is left with a 3-dimensional 
Schrodinger equation depending on the inter-particle distances only. We use the perimetric 
coordinates 

x = n + r 2 - r 3 , 

y = r 1 -r 2 +r 3 , (6) 
z = —r\ +r 2 + r 3 , 

and express the Schrodinger equation as a generalized eigenvalue problem for the energy E: 

A\V) =E B (7) 

where A and B operators are polynomials in the x,y,z,d x ,d y ,d z operators. The basis func- 
tions used in the calculation are direct products of Laguerre polynomials and exponentials 
along each perimetric coordinate, whose properties are discussed in details in 0. The matri- 
ces representing A and B in such a basis set are real symmetric sparse banded matrices, where 
all matrix elements are known in analytic form and involve only simple algebraic expressions. 
The generalized eigenvalue problem is then solved using the Lanczos algorithm in order to 
produce few eigenvalues, among several thousand, in the interesting energy range. 

Whereas an accurate computation of the ground v = 0, L = level of the 2ptr„ state 
requires only a moderately large basis set (about 20,000), the computation of the first excited 
state is much more difficult and requires at least a basis size of 150,000 and a careful choice of 
the variational parameters of the basis. We used basis sizes up to 450,000 to confirm that the 
results discussed here are fully converged. For the first excited vibrational level we find a total 
energy E = -0.499727840801511 a.u. and a dissociation energy E = -0.499727839716466 a.u. 
This gives a binding energy B = 1.085045 x 10~ 9 a.u. with an uncertainty of the order of 
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r (Bohr radii) 

Fig. 2 - Wavefunctions (not normalized) of the ground (dashed line) and excited (solid thick line) 
levels of the 2pa u state of the molecular ion. The last one is compared with the corresponding 
p-H zero energy scattering wavefunction (thin line) . The existence of an excited level with very small 
binding energy is predicted by our calculations. Its wavefunction extends very far in the internuclear 
distance r, with a maximum probability density around 100 Bohr radii. It is responsible for a huge 
scattering length of 750 a.u. Note the use of a logarithmic scale on r. 

We have also calculated the wavefunctions. These are full three-body wavefunctions and 
are thus not easily plotted. However, they take significant values only for rather large in- 
ternuclear distances: the electronic wavefunction is thus essentially the ground state of the 
hydrogen atom attached to one of the protons, independently of the internuclear distance. We 
checked that this simple property is almost exactly obeyed by the three-body wavefunction. 
Once this trivial part of the wavefunction is factored out, one is left with a wavefunction 
depending only on the internuclear distance. It is plotted in fig. for the lowest two levels. 
Because of the very large size of the excited state, we chose to plot the wavefunction using a 
logarithmic scale for the internuclear distance r. The ground level is a nodeless wavefunction 
centered around r = 15 a.u., while the excited level extends much further; it has a maximum 
at r k 100 a.u. and still significant values at r ~ 1000 a.u. There is an inflexion point at 
r ~ 215 a.u., located at the outer turning point of the Born-Oppcnhcimer (or Landau) po- 
tential, where the 2pa u potential equals the binding energy. In the same figure, we also plot 
the zero-energy wavefunction got from the Faddeev approach. At small r, it is remarkably 
similar to the wavefunction of the excited level, which is not surprising considering the very 
small energy difference. At large distance, the zero-energy wavefunction diverges linearly and 
has a zero at r » a, the scattering length value 750 a.u. 

Results presented above have been obtained with the mass ratio m p /m e = 1836.152701, 
as recommended in |16|. Using the more recent 1998 CODATA value fjT] would not change 
anything significantly. Our method makes possible to compute the energy levels for any mass 
ratio of the particles. When the mass ratio is decreased, the binding energy also decreases 
until a critical value beyond which the excited level disappears and one is left with a single 
bound state. We estimate this critial mass ratio to be around 1781. This is rather close 
to the actual value which explains why the state is so weakly bound. The closeness of the 
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critical mass ratio also explains why the Landau potential discussed above - although it is 
fairly accurate - gives a wrong binding energy. 

Calculations presented here were performed using a fully non relativistic dynamics with 
just Coulomb pair-wise interactions taken into account. In view of the extreme sensitivity of 
these results, it is necessary to quantify the possible relativistic effects. 

As the dynamics is governed by a shallow potential well at distances around 10 a.u., 
where the nuclear motion is very slow, relativistic corrections must be considered only for 
the electronic motion. As the internuclear distance is much larger than the typical electron- 
nucleus distance, relativistic corrections will be essentially identical for the weakly bound state 
and the dissociation limit, leaving the binding energy only very weakly modified. The first 
order relativistic and radiative corrections for the states have been obtained in |18j and 
are discussed with some detail in jS] . The results summarizing the relative corrections to the 
binding energy (AB / B) for the L = lsc g (filled circles) and 2pa u (filled square) states are 
given in fig. O For lscg states they are smaller than 10 -3 and vary smoothly over five decades 
of binding energy. For the v = 2pa u level they are one order of magnitude smaller. In 
absolute value they are AB ~ 10~ 9 a.u. and as they scale as 1/r 4 [18] . they should be at least 
one or two orders of magnitude smaller for v — 1. Thus, despite the smallness of its binding 
energy, relativistic effects preserve the bound character of the new state we have presented. 
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Fig. 3 - Relativistic corrections to the binding energies for the L — ls<r fl (circles) and the ground 
2pa u (square) state of the H^" molecular ion. 



Another source of relativistic corrections, not included in Refs. |18ll3). is the modification 
of the 1/r 4 polarization potential at very large distances due to retardation effects (Casimir- 
Polder effect). Contrary to the dipole-dipole case (Van der Waals forces), where the 1/r 6 
long-range behaviour is changed into 1/r 7 , in the charge-dipole case we are considering the 
usual 1/r 4 term is modified by adding a 1 /r 5 contribution [19] in the form: 



V(r) 



ad 
~2r 4 



1 - 



11a m e 1 
2tt m p r 



(8) 



where a is the fine structure constant. At r 
negligible. 



100 a.u., the correction turns out to be 
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One should be careful with the spin-orbit and spin-spin interactions. As all states consid- 
ered here have zero total angular momentum L, the spin-orbit coupling vanishes and all states 
correspond to J = 1/2 (where J is the sum of the total angular momentum L and the electronic 
spin S. The spin-spin interaction is more tricky. Indeed, we are interested in triplet states with 
total two-proton spin S pp = 1, so that the total angular momentum can be either F — 1/2 
or F — 3/2. This hyperfme structure should be very close to the F = / F = 1 hyperfme 
structure of the hydrogen atom. As the latter is much larger than the binding energy we have 
calculated for the non-relativistic problem, it is likely that the v = 1, L = 0, J = 1/2, F — 3/2 
level lies above the dissociation limit of the F = 1/2 series. The dissociation rate induced by 
the hyperfme coupling is however most probably very low. 

A direct measurement of the p-H cross section at very low energy seems unlikely at present. 
One can however access the low energy p-H continuum in the final state of the photodis- 
sociation cross section. The excited vibrational 2p<x u (v = 1,L = 0) level predicted here is 
radiatively coupled to the lsa g (v — 19, L = 1) level. The electric dipole transition between 
those two levels should be observable in the 6 GHz range using an experiment similar to the 
one used to detect the (v — 0, L = 0) — > (v = 19, L = 1) transition [20112 1| . An experimental 
confirmation of our results would be very interesting. 

* * * 
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